SPECIFICATION 



TO ALL WHOM IT MAY CONCERN: 

BE IT KNOWN that we, MASAHIRO ISHIDA a subject of Japan and 
residing at Nerima-ku, Tokyo, Japan, TAKAHIRO YAMAGUCHI a subject 
of Japan and residing at Nerima-ku, Tokyo, Japan and MANI SOMA a 
subject of U.S.A. residing at Seattle, WA U.S.A. have invented certain new 
and useful improvements in 

"APPARATUS FOR AND METHOD OF MEASURING JITTER" 

and we do hereby declare that the following is a full, clear and exact 
description of the same; reference being had to the accompanying drawings 
and the numerals of reference marked thereon, which form a part of this 
specification. 



-1- 

FOR AND METHOD OF MEASURING JITTER 

BACKGROUND OF THE INVENTION 

The present invention particularly relates to an apparatus for and a 
5 method of measuring period jitter that are applied to a measurement of jitter 
of, for example, a microprocessor clock or another periodic signal. 

A time interval analyzer and/or an oscilloscope have conventionally 
been used for the measurement of period jitter. The method of these 
apparatus is called Zero-crossing Method, in which, as shown in Fig. 1, a 

10 clock signal (a signal under measurement) x(t) from, for example, a PLL 

(Phase-Locked Loop) under test 1 1 is supplied to a time interval analyzer 12. 
Regarding a signal under measurement x(t), a next rising edge following one 
rising edge fluctuates against the preceding rising edge as indicated by dotted 
lines. That is, a time interval T p between the two rising edges, namely a 

1 5 period fluctuates. In the Zero-crossing Method, a time interval between 
zero-crossings (period) of the signal under measurement is measured, and a 
fluctuation of period is measured by a histogram analysis. The histogram is 
displayed, for example as shown in Fig. 2, by the time interval analyzer 12. 
A time interval analyzer is described in, for example, "Phase Digitizing 

20 Sharpens Timing Measurements" by D. Chu, IEEE Spectrum, pp. 28-32, 
1988, and "A Method of Serial Data Jitter Analysis Using One-Shot Time 
Interval Measurements" by J. Wilstrup, Proceedings of IEEE International 
Test Conference, pp. 819-823, 1998. 

Tektronix, Inc. and LeCroy Co. have recently been providing digital 

25 oscilloscopes each being able to measure a jitter using an interpolation 

method. In this jitter measurement method using the interpolation method 
(interpolation-based jitter measurement method), an interval between data 
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having signal values close to a zero-crossing out of measured data of a 
sampled signal under measurement is interpolated to estimate a timing of 
zero-crossing. That is, in order to estimate a fluctuation of period, a time 
interval between zero-crossings (period) is estimated using a data 
5 interpolation with a small error. 

That is, as shown in Fig. 3, a signal under measurement x(t) from the 
PLL under test 1 1 is inputted to a digital oscilloscope 14. In the digital 
oscilloscope 14, as shown in Fig. 4, the inputted signal under measurement 
x(t) is converted into a digital data sequence by an analog-to-digital converter 

10 15. A data-interpolation is applied to an interval between data having signal 
values close to a zero-crossing in the digital data sequence by an interpolator 
16. With respect to the data-interpolated digital data sequence, a time 
interval between zero-crossings is measured by a period estimator 17. A 
histogram of the measured values is displayed by a histogram estimator 18. 

15 In addition, a root-mean-square value and a peak-to-peak value of fluctuations 
of the measured time intervals are obtained by an RMS & Peak-to-Peak 
Detector 19. For example, in the case in which a signal under measurement 
x(t) has a waveform shown in Fig. 5A, its period jitters are measured as 
shown in Fig. 5B. 

20 On the other hand, inventors of the present invention have proposed a 

method of measuring a jitter as described below in an article entitled 
"Extraction of Peak-to-Peak and RMS Sinusoidal Jitter Using an Analytic 
Signal Method" by TJ. Yamaguchi, M. Soma, M. Ishida, and T. Ohmi, 
Proceedings of 1 8th IEEE VLSI Test Symposium, pp.395-402, 2000. That 

25 is, as shown in Fig. 6, an analog clock waveform from a PLL (Phase locked 
loop) circuit under test 1 1 is converted into a digital clock signal x c (t) by an 
analog-to-digital converter 22, and the digital clock signal x c (t) is supplied to 



a Hilbert pair generator 23 acting as an analytic signal transforming part, 
where the digital clock signal x c (t) is transformed into an analytic signal Zc(t). 

In the Hilbert pair generator 23, signal components around a 
fundamental frequency of the clock signal under measurement x c (t) are 
extracted by a band-pass filter (not shown), and the extracted signal 
components are inputted to a Hilbert transformer 25, where the signal 
components are Hilbert-transformed. An analytic signal z c (t) having the 
signal component not Hilbert-transformed and the signal component 
Hilbert-transformed as a real part and an imaginary part of a complex number, 
respectively is obtained from the Hilbert pair generator 23. 

An instantaneous phase 0(t) of the clock signal x c (t) is estimated by 
an instantaneous phase estimator 26 from the analytic signal z c (t). A linear 
phase is removed by a linear phase remover 27 from this instantaneous phase 
0(t), and an instantaneous phase noise waveform A<f>(t) is obtained. 

This instantaneous phase noise A<j)(t) is sampled by a zero-crossing 
sampler 28, and the sampled instantaneous phase noise is inputted to a 
peak-to-peak detector 32 as a timing jitter sequence A$[n], where a difference 
between the maximum peak value max (A(j)[k]) and the minimum peak value 
min (A(|>[k]) of the A<()[n] (=A(j>(nT)) is calculated to obtain a peak value 
(peak-to-peak value) A(j) pp of timing jitter as follows. 

max (A(j)[k])- min (A<j)[k]) 

k k - 4 

In addition, the timing jitter sequence A<))[n] is also inputted to a 
root-mean-square detector 33, where a root-mean-square (RMS) value of the 
timing jitter sequence A<|)[n] is calculated using following equation to obtain a 
root-mean-square value AcJjrms of timing jitter. 



1 N ~ l 

A(j>RMs= — S A ^ 2 W 

This method is referred to as the A<j) method, since a peak value of 
timing jitter (peak-to-peak value) and a root-mean-square value of timing 
5 jitter are obtained as described above from the instantaneous phase noise 
waveform A<j)(t). 

According to the A(j> method, a timing jitter can be measured at high 
speed with relatively high accuracy. 

Since the jitter measurement method in the time interval analyzer 

10 system includes an intermediate dead-time during which no measurement can 
be performed after each one of periodic measurements is performed, there is a 
problem that it takes a long time to acquire a number of data that are required 
for a histogram analysis. 

In addition, in a jitter measurement method in which a wide-band 

15 oscilloscope and an interpolation method are combined, there is a problem 
that a jitter is overestimated (overestimation). That is, the measured jitter 
values in this method are not compatible with the values measured by the 
zero-crossing method. For example, a jitter measurement result by the time 
interval analyzer and a jitter measurement result by the interpolation method 

20 for a 400 MHz clock signal are shown in Figs 7 A and 7B, respectively so that 
those measurement results can be compared with each other. According to 
those figures, the measured value by the time interval analyzer is 7.72 ps 
(RMS), while the measured value by the interpolation method is 8.47 ps 
(RMS), and this value is overestimated. In addition, the jitter measurement 

25 method by the interpolation method requires a long measurement time. 

In addition, in the conventional A§ method, a zero-crossing point of a 
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signal under measurement is approximated by a sampling point closest to 
each zero-crossing point (this point is referred to as adaptive zero-crossing 
points approximation). Therefore, there is a possibility that an amount of 
error in jitter measurement is increased when an over-sampling rate in the 
5 signal digitization process is small. 

It is an object of the present invention to provide a jitter measurement 
apparatus and its method that can estimate a jitter value compatible with a 
jitter value measured by the conventional time interval analyzer method. 

10 SUMMARY OF THE INVENTION 

The jitter measurement apparatus according to the present invention 
comprises: phase error estimation means for estimating sampling points close 
to corresponding zero-crossing points of the signal under measurement as 
approximated zero-crossing points, for obtaining phase errors between the 

15 approximated zero-crossing points and the corresponding zero-crossing points 
of the signal under measurement, and for outputting a zero crossing phase 
error data sequence; and a period jitter estimation means for obtaining a 
period jitter sequence of the signal under measurement from the phase error 
data sequence. 

20 It is desirable that the jitter measurement apparatus according to the 

present invention further comprises: cycle-to-cycle period jitter estimation 
means to which the period jitter sequence is inputted for calculating its 
difference sequence, and for outputting a cycle-to-cycle period jitter 
sequence. 

25 It is desirable that the jitter measurement apparatus according to the 

present invention further comprises: band-pass filter means for passing 
therethrough predetermined frequency components of the signal under 
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measurement, and for supplying the predetermined frequency components to 
the phase error estimation means. 

In addition, it is desirable that the jitter measurement apparatus 
according to the present invention further comprises: jitter detection means to 
5 which a jitter sequence is inputted for obtaining a jitter value of the signal 
under measurement from the jitter sequence. 

It is desirable that the phase error estimation means comprises: 
zero-crossing sampling means to which the signal that has passed through the 
band-pass filter means is inputted for sampling only waveform data around 
10 zero-crossing timings of the band-limited signal, and for outputting an 

approximated zero-crossing data sequence; and phase error calculation means 
to which the approximated zero-crossing data sequence is inputted for 
calculating phase errors between the approximated zero-crossing points and 
the corresponding zero-crossing points of the signal under measurement from 
15 the approximated zero-crossing data, and for outputting the zero-crossing 
phase error data sequence. 

It is desirable that the phase error estimation means further comprises: 
analytic signal transformation means to which the signal that has passed 
through the band-pass filter means is inputted for transforming the 
20 band-limited signal into a complex analytic signal, and for supplying the 
complex analytic signal to the zero-crossing sampling means. 

It is desirable that the period jitter estimation means comprises: 
instantaneous period estimation means to which the zero-crossing phase error 
data sequence is inputted for obtaining an instantaneous period sequence of 
25 the signal under measurement from sampling intervals of the approximated 
zero-crossing data sequence and the zero-crossing phase error data; 
fundamental period subtraction means for obtaining differences between the 
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instantaneous periods and a fundamental period of the signal under 
measurement, and for outputting an uncorrected period jitter sequence; and 
period jitter sequence correction means for multiplying the uncorrected period 
jitter sequence by a ratio of the fundamental period of the signal under 
5 measurement and the zero-crossing sampling intervals to perform a correction 
of the uncorrected period jitter sequence. 

It is desirable that the period jitter estimation means comprises: 
instantaneous angular frequency estimation means to which the zero-crossing 
phase error data sequence is inputted for obtaining an instantaneous angular 

10 frequency sequence of the signal under measurement from the sampling 
intervals of the approximated zero-crossing data sequence and the 
zero-crossing phase error data; and period jitter calculation means for 
obtaining the period jitter sequence from the instantaneous angular frequency 
sequence and the fundamental period. 

15 It is desirable that the band-pass filter means comprises: time domain 

to frequency domain transforming means for transforming the signal under 
measurement into a signal in frequency domain; bandwidth limitation means 
for extracting only components around the fundamental frequency of the 
signal under measurement from an output of the time domain to frequency 

20 domain transformation means; and frequency domain to time domain 

transformation means for inverse-transforming an output of the bandwidth 
limitation means into a signal in time domain. 

The jitter detection means comprises one or a plurality of 
peak-to-peak detection means for obtaining a difference between the 

25 maximum value and the minimum value of a supplied jitter sequence, RMS 
detection means for obtaining a root-mean-square value (RMS value) of a 
supplied jitter sequence, and histogram estimation means for obtaining a 
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histogram of a supplied jitter sequence. 

The jitter measurement method according to the present invention 
comprises: a step of estimating sampling points close to corresponding 
zero-crossing points of the signal under measurement as approximated 
5 zero-crossing points, obtaining phase errors between the approximated 

zero-crossing points and the corresponding zero-crossing points of the signal 
under measurement, and estimating a zero-crossing phase error data sequence; 
and a step of obtaining a period jitter sequence of the signal under 
measurement from the phase error data sequence. 
1 0 In addition, it is desirable that the jitter measurement method 

according to the present invention further comprises a step of calculating a 
difference sequence of the period jitter sequence to obtain a cycle-to-cycle 
period jitter sequence of the signal under measurement. 

It is desirable that the jitter measurement method according to the 
1 5 present invention further comprises a step of passing predetermined frequency 
components of the signal under measurement to band-limit the signal under 
measurement, and moving to the step of estimating the zero-crossing phase 
error data sequence. 

It is desirable that the jitter measurement method according to the 
20 present invention further comprises a step of obtaining a jitter value of the 
signal under measurement from a jitter sequence. 

It is desirable that the step of estimating the phase error data sequence 
comprises: a step of sampling only waveform data close to zero-crossing 
timings of the band-limited signal, and estimating an approximated 
25 zero-crossing data sequence; and a step of calculating phase errors between 
the approximated zero-crossing points and the corresponding zero-crossing 
points of the signal under measurement from the approximated zero-crossing 
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data sequence, and obtaining the zero-crossing phase error data sequence. 

It is desirable that the step of estimating the phase error data sequence 
comprises: a step of transforming the band-limited signal into a complex 
analytic signal; a step of sampling the analytic signal waveform data close to 
5 zero-crossing timings of the transformed analytic signal to estimate an 

approximated zero-crossing complex data sequence; and a step of calculating 
phase errors between the approximated zero-crossing points and the 
corresponding zero-crossing points of the signal under measurement from the 
approximated zero-crossing complex data sequence, and obtaining the 
1 0 zero-crossing phase error data sequence. 

It is desirable that the step of obtaining the period jitter sequence 
comprises: a step of obtaining an instantaneous period sequence of the signal 
under measurement from the sampling intervals of the approximated 
zero-crossing data sequence and the zero-crossing phase error data; a step of 
15 calculating differences between the instantaneous periods and the 

fundamental period of the signal under measurement, and obtaining an 
uncorrected period jitter sequence; and a step of multiplying the uncorrected 
period jitter sequence by a ratio of the fundamental period of the signal under 
measurement and the corresponding zero-crossing sampling intervals to 
20 correct the uncorrected period jitter sequence values. 

It is desirable that the step of obtaining the period jitter sequence 
comprises: a step of obtaining an instantaneous angular frequency sequence of 
the signal under measurement from the sampling intervals of the 
approximated zero-crossing data sequence and the zero-crossing phase error 
25 data; and a step of obtaining a period jitter sequence from the instantaneous 
angular frequency sequence and the fundamental period. 

It is desirable that the step of band-limiting the signal under 
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measurement comprises: a step of storing the signal under measurement in a 
buffer memory; a step of extracting the signal in the sequential order from the 
buffer memory such that each partial signal being extracted is partially 
overlapped with the partial signal extracted just before; a step of multiplying 
5 each extracted partial signal by a window function; a step of transforming 
each partial signal multiplied by the window function into a spectrum signal 
in frequency domain; a step of extracting only components around the 
fundamental frequency of the signal under measurement from the spectrum 
signal; a step of inverse-transforming the spectrum signal having only the 

10 fundamental frequency component into a signal in time domain; and a step of 
multiplying the inverse-transformed signal in time domain by an inverse 
number of the window function to obtain the band-limited signal. 

A principle of the present invention will be described below. In this 
description, a clock signal is used as a signal under measurement. However, 

15 the present invention can be applied to any periodic signal. 



JITTER MEASUREMENT METHOD 

A jitter-free clock signal is a square wave having a fundamental 
frequency f 0 . This signal can be decomposed by Fourier analysis to 
20 harmonics having frequencies f 0 , 3f 0 , 5f 0 , Since a jitter corresponds to a 
fluctuation of a fundamental frequency of a signal under measurement, only 
signal components around the fundamental frequency are handled in the jitter 
analysis. 

A fundamental sinusoidal wave component x(t) of a jittery clock 
25 signal (clock signal under measurement) can be expressed, assuming that its 
amplitude is A and its fundamental period is To, by the following equation 

(i). 
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x(t)=Acos((j)(t))=Acos((27u/To)t+e-A(|)(t)) 



(1) 



Here, (|)(t) is an instantaneous phase of the signal under measurement, and can 
be expressed by a sum of a linear instantaneous phase component 2nt/T 0 
5 containing the fundamental period T 0 , an initial phase angle component 0, and 
an instantaneous phase noise component A<|)(t). 

An analytic signal z(t) of a real signal x(t) is defined as a complex 
signal expressed by the following equation (2). 



In this case, j represents an imaginary unit, and an imaginary part x (t) of the 
complex signal z(t) is a Hilbert transform of a real part x(t). 

On the other hand, the Hilbert transform of a time waveform x(t) is 
1 5 defined by the following equation (3). 



In this case, x (t) is a convolution of the function x(t) and (1/jcf). That is, the 
20 Hilbert transform is equivalent to an output at the time when the x(t) is passed 
through an entire-band-pass filter. However, in the output x (t) in this case, 
the magnitude of its spectrum components has not been changed, but its phase 
has been shifted by n/2. Analytic signal and Hilbert transform are described 
in, for example, "Probability, Random Variables, and Stochastic Processes" 
25 by A. Papoulis, 2nd edition, McGraw-Hill Book Company, 1984. 



10 



z(t) = x(t)+jx(0 



(2) 




(3) 



Therefore, an analytic signal of a fundamental cosine wave x(t) of a 



1 
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signal under measurement is given by the following equation (4). 

Z(t)=x(t)+jx(t) 

= Acos((27u/T 0 )t+e-A())(t)) 
+jAsin((27c/T 0 )t+e-A(j)(t)) (4) 

An instantaneous frequency of the analytic signal z(t) is given, as is well 
known, by the following equation (5). 

1 ©(t) 1 x(t)H'[x(t)]-x'(t)H[x(t)] 



T 0 +J 2n 2% x 2 (t) + H 2 [x(t)] 
10 [Hz] (5) 



1 fl-^_A<j>'(t)" 
{ 2n , 



T 0 



That is, a period jitter J can be obtained by the following equation (7) from an 
equation (6). 



15 T 0 +J(t)*T 0 fl + ^A<|>'(t)| [sec] (6) 

^ 2U J 

T 2 

J(t)*-2-Af(t) [sec] (7) 

2tz 



On the other hand, an instantaneous angular frequency co(t) is given 
by the following equation (8). 

20 

co(t)- ((2Tc/T 0 )t+e-A(j)(t)) , -(27i/To)-A<l) t (t) [rad/sec] (8) 



I 
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Therefore, from the equations (7) and (8), a period jitter J(t) can be expressed 
by the following equation (9) using an instantaneous angular frequency co(t). 

5 J(tMT 0 2 /2tt^ [sec] (9) 

As shown in Fig. (8), when adjacent rising zero-crossing points of a 
signal under measurement x(t) are approximated by corresponding sampling 
points t[k] and t[k+l] that are closest to the corresponding zero-crossing 
10 points t k and t k +i, an instantaneous angular frequency oo(t) between the two 

adjacent approximated zero-crossing points [t[k], t[k+l]] is given by a ratio of 
a phase angle shifting during a time between t[k] and t[k+l] and a length of 
the time interval, assuming that a phase error between the approximated 
zero-crossing point t[k] and the real zero-crossing point t k is 8[k]. 

15 

co[k]-{(27i+(6[k+l]-5[k]))/T k?k+1 [rad/sec] (10) 

In this case, the T^k+i represents the length of the time interval between the 
two adjacent approximated zero-crossing points t[k] and t[k+l], and can be 
20 obtained, as shown in Fig. 8, by obtaining a timing sequence t[k] of the 
approximated zero-crossing points and then calculating differences of the 
timing sequence t[k], 

T u+ i=t[k+l]-t[k] [sec] (11) 

25 

In addition, a phase error §[k] can be expressed by a difference 
between a phase angle cj>(t[k]) at an approximated zero-crossing point t[k] and 



I 
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a phase angle <|>(tk) at a zero-crossing point t k . 

5[k]= Kt[k]> ())(t k ) [rad] (12) 

5 Here, a band-limited signal under measurement x(t) is a fundamental cosine 
wave having an amplitude A, as shown in the equation (1), and a phase 
angle <|)(t) and an amplitude value x(t) are in the following relationship. 

x(t)=Acos((()(t))=Asin((t)(t)+(7r/2)) [rad] (13) 

10 

For this reason, a phase angle (j)(t) of a signal under measurement at a time 
point t can be expressed by the following equation (14) from the amplitude 
value x(t) at the time point t and the amplitude A using an arcsine function. 

15 ^(t^sin^CxCtyAK^) t rad ] ( 14 ) 

However, since an arcsine function can be expressed using a phase value in 
the range between -te/2 and +7i/2, the equation (14) is valid in a range between 
-7i and 0 of the phase angle of the cosine wave, i.e., in a rising region of the 
20 cosine wave. Therefore, if an amplitude value at an approximated 
zero-crossing point t[k] is assumed as expressed by an equation (15) 

x Z x[k]=x(t[k]) (15) 

25 regarding a rising zero-crossing point of a fundamental cosine wave x(t), a 
phase angle §(t[k]) at t[k] is given by the following equation (16) from the 
equation (14). 
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(t)(t[k])=sm I (xzx[k]/A)-(7u/2) [rad] (16) 

On the other hand, since an amplitude value at a real rising zero-crossing 
5 point t k is 0, a phase angle <j)(t k ) at t k is given by the following equation (17). 

(J)(t k )=sm 1 (0/A)-(7i/2)=-(7i/2) [rad] ( 1 7) 

Therefore, a phase error 5[k] between the rising zero-crossing point t k and its 
10 approximated zero-crossing point t[k] can be obtained by the following 
equation (18) from the equations (12), (16) and (17). 

5[k]=arcsin(x zx [k]/A) [rad] (1 8) 

15 In addition, in the range between 0 and n of the phase angle of the cosine 

wave, i.e., in a falling region, a phase angle <|>(t) at a time point t is expressed 
by the following equation (19) from the amplitude value x(t) at the time point 
t and the amplitude A using an arcsine function. 

20 (f>(t)= -(sin ! (x(t)/A)-(7c/2)) [rad] (19) 

Therefore, regarding a falling zero-crossing point of a fundamental cosine 
wave x(t), a phase error 5[k] between the approximated zero-crossing point 
and the real zero-crossing point can be obtained by the following equation 
25 (20). 



8 [k]—arcsin(x zx [k] / A) 



[rad] (20) 
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In addition, a phase error 5[k] may be obtained by transforming a 
signal under measurement x(t) to an analytic signal z(t), 

5 z(t)=x(t)+jx(t) (21) 

then by obtaining, from an analytic signal data z(t[k]) at the approximated 
zero-crossing point, 

1 0 z(t[k])=x(t[k])+j x (t[k]) s x zx [k]+j x zx [k] (22) 

an instantaneous phase (|)[k] at the approximated zero-crossing point, 

(Kk^an^xzxM/xzxM] [rad] (23) 

15 

and then by calculating a difference from the instantaneous phase <|>[k] at the 
approximated zero-crossing point. 

That is, regarding the rising zero-crossing point, a difference 5[k] between the 
instantaneous phase (|)[k] at the approximated zero-crossing point and the 
20 instantaneous phase -%/2 of the rising zero-crossing point is obtained. 

5[k]=(|)[k]+(7i/2) [rad] (24) 

Regarding the falling zero-crossing point, a difference 8[k] between the 
25 instantaneous phase <()[k] at the approximated zero-crossing point and the 
instantaneous phase n/2 of the rising zero-crossing point is obtained. 
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5[k]=nj>[k]-(7c/2) 



[rad] 



(25) 



By the processes described above, a period jitter J[k] can be obtained 
by the following equation (26) from the equations (9) and (10) using the 
5 zero-crossing phase error data sequence 8[k]. 



J[k] = 



L k,k+1 



§[k + l]-5[k] 

271 



-T ft 



T, 



k,k+l 



[sec] (26) 



In this case, the terms in the parentheses() in right side represent, as shown in 
10 Fig. 8, a difference between an instantaneous period T[k] 



T[k] = T k?k+1 



5[k + l]-5[k] 
To" 



[sec] 



(27) 



of the signal under measurement and the fundamental period T 0 , i.e., a period 
1 5 j itter (uncorrected period jitter). In addition, the term of (T 0 /T k , k+ i) is a 
correction term to the period jitter obtained from the approximated 
zero-crossing point. A period jitter can be obtained with higher accuracy by 
this correction. 

As shown in Fig. 9, the correction term can reduce an error between 
20 an estimated value and a theoretical value of jitter with respect to an RMS 
value and a peak-to-peak value of period jitter. That is, when the correction 
term is not used, there are errors of +0.53% and +15.9% to the RMS value 
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and the peak-to-peak value, respectively, but in the method of the present 
invention using the correction term, those errors are remarkably decreased to 
-0.004% and +0.04%, respectively. Further, the data in Fig. 9 are obtained 
from a computer simulation with respect to a signal having a sinusoidal jitter. 
5 In addition, an experimental result using a real waveform is shown in Fig. 
1 OA. The numbers of samplings per period are shown in the lateral axis, and 
RMS values of period jitter are shown in the longitudinal axis. In this case, 
the RMS value is substantially constant regardless of the number of samplings, 
and similar results are obtained in the case where the correction term is used 

1 0 and in the case where the correction term is not used. On the other hand, it 
is understood from an experimental result that the method not using the 
correction term overestimates peak-to-peak values of period jitter as the 
over-sampling rate decreases, but peak-to-peak values of period jitter can 
accurately be obtained by the method of the present invention using the 

1 5 correction term. Particularly, when the over-sampling rate is small, the 

effect of the present invent is large. For example, in the example of Fig. 10B, 
according to the present invention, when the number of samplings per period 
is 8 points (over-sampling rate is 4), the error can be corrected by 
approximately 8%, and when the number of samplings per period is 3 points 

20 (over-sampling rate is 1.5), the error can be corrected by approximately 1 8%. 

In the jitter measurement method according to the present invention, 
in order to extract a fundamental frequency component of the clock signal 
under measurement x(t), the clock signal under measurement x(t) is 
band-limited first using a band-pass filter. The band-limited clock signal 

25 x BP (t) is sampled at the timings (approximated zero-crossing points) closest to 
the corresponding zero-crossing points of the band-limited clock signal x BP (t), 
then a phase error sequence 8[n] between the approximated zero-crossing 
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points and the corresponding zero-crossing points of the signal under 
measurement is obtained from this zero-crossing sample sequence x zx [n]. 
Next, a period jitter sequence of the signal under measurement is estimated 
using the formula in parentheses () of the equation (26) from the phase error 
5 sequence 5[n]. In this case, the period for obtaining a period jitter may be 
m periods (m=0.5, 1, 2, 3, • • •)• F° r example, by taking the m=0.5 period, a 
fluctuation of a time interval between a rising (or falling) zero-crossing point 
and a next falling (or rising) zero-crossing point may be obtained. Also, by 
taking the m=2 periods, a fluctuation of a time interval between a rising (or 

1 0 falling) zero-crossing point and a rising (or falling) zero-crossing point after 
two rising (or falling) points from the first rising (or falling) zero-crossing 
point may be obtained. Then the period jitter sequence data are multiplied 
by the correction term (T 0 /T Kk+l ) to correct the values of the period jitter 
sequence. Fig. 1 1 shows a period jitter waveform J[n] of the signal under 

1 5 measurement x(t) obtained using m=l period. In addition, by calculating a 
root-mean-square value and a difference value between the maximum value 
and the minimum value of the period jitter data measured in this manner using 
the following equations (28) and (29) respectively, an RMS Jrms value and a 
peak-to-peak J P p value of the period jitter can be obtained. 

20 

j - = JM j2[k] [sec] (28) 

J PP =max(J[k])-min(J[k]) [sec] (29) 

k k 



25 
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In this case, M is the number of samples of the measured period jitter data. 
Fig. 12 shows a histogram (Fig. 12A) of the period jitter measured by the 
jitter measurement method according to the present invention and a histogram 
(Fig. 12B) measured by the conventional time interval analyzer. Thus, the 
5 histogram measured by the method of the present invention can be 

comparable with the histogram of the conventional time interval analyzer. 
Both histograms indicate the similar waveforms to each other. In addition, 
Fig. 13 shows an RMS value and a peak-to-peak value of the period jitter 
measured by the jitter measurement method according to the present invention 

10 as well as the corresponding values measured by the time interval analyzer. 
In this case, the peak-to-peak value J PP of the observed period jitter is 
substantially proportional to a square root of logarithm of the number of 
events (the number of zero-crossings). For example, in the case of 
approximately 5000 events, J PP =45 ps is a correct value. A J PP error in Fig. 

15 13 is shown assuming that 45 ps is the correct value. As shown in Figs. 
12 A, 12B, and 13, the jitter measurement method according to the present 
invention can provide measured jitter values comparable with the 
conventional measurement method. 

Furthermore, the jitter measurement method according to the present 

20 invention can measure a cycle-to-cycle period jitter and a period jitter 

simultaneously. A cycle-to-cycle period jitter J C c is a period fluctuation 
between adjacent cycles, and is expressed by the following equation (30). 



25 



JccM=T[k+l]-T[k] 

=(T 0 +J[k+1])-(T 0 +J[k]) 
=J[k+l]-J[k] [sec] 



(30) 
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Therefore, by obtaining differences between the period jitter data measured as 
described above, and then by calculating their root-mean-square and a 
difference between the maximum value and the minimum value using the 
following equations (31) and (32) respectively, an RMS value J C c,rms and a 
5 peak-to-peak value J C c,pp of cycle-to-cycle period jitter can be obtained. 



Jcc 5 RMS=^ZJccM [sec] (31) 

J CC)PP =max(J cc [k])-rnin(J cc [k]) [sec] (32) 

10 

In this case, L is the number of samples of the measured cycle-to-cycle period 
jitter data. Fig. 14 and Fig. 15 show waveforms of a cycle-to-cycle period 
jitter J cc [n] and its histogram obtained by the jitter measurement method 
according to the present invention, respectively. 

1 5 The j itter measurement method according to the present invention can 

be applied, as described above, not only to the estimation of jitter of a clock 
signal but also to the estimation of jitter of the other signals. 

In the jitter measurement method according to the present invention, a 
band-pass filtering process may be applied, after an analog signal under 

20 measurement is digitized, to the digitized signal under measurement, or the 
analog signal waveform may be digitized, after the band-pass filtering process 
is applied to the analog signal under measurement. The band-pass filtering 
process may be performed by either one of an analog filter, a digital filter and 
a software using Fourier transform. In addition, it is desirable that a high- 

25 speed AD converter, a high-speed digitizer or a high-speed digital sampling 
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oscilloscope is used for the digitization of the analog signal. 

In addition, in the jitter measurement method according to the present 
invention, a period jitter can also be measured with high accuracy by 
removing amplitude modulation (AM) components of the signal under 
5 measurement using waveform clipping means in the state that phase 
modulation (PM) components corresponding to jitter are retained. 

BAND-PASS FILTER 

A band-limitation of a digitized digital signal can be realized by a 

1 0 digital filter, and can also be performed by Fourier transform. Next, a 
band-pass filtering using FFT (Fast Fourier Transform) will be explained 
below. FFT is a method of transforming at high speed a signal waveform in 
time domain to a signal in frequency domain. 

First, a digitized signal under measurement x(t) shown in Fig. 16 is 

1 5 transformed to a signal X(f) in frequency domain using FFT. Fig. 1 7 shows 
a power spectrum of the transformed signal X(f). Then data around a 
fundamental frequency (400 MHz in this example) of the signal X(f) are 
retained and the other data are made zero. Fig. 18 shows a signal X BP (f) of 
the band-limited signal in frequency domain. Finally, a band-limited signal 

20 waveform x BP (t) in time domain can be obtained by applying inverse FFT to 
the band-limited signal X BP (f) in frequency domain. Fig. 19 shows the 
band-limited signal waveform x B p(t) in time domain. 

APPROXIMATED ZERO-CROSSING POINT DETECTION METHOD 
25 Next, an approximated zero-crossing point detection method will be 

described. First, the maximum value of a real part x(t) of an analytic signal 
of an inputted signal under measurement is defined as 100% level, and the 
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minimum value is defined as 0% level to calculate 50% level signal value V 
(50%) as a zero crossing level. Then, a difference between a sample value 
x(j-l) and the 50% level V (50%) of the x(t) and a difference between its 
adjacent sample value x(j) and the 50 % level V(50%) of the x(t), i.e., 
5 (x(j-l)-V(50%)) and ( x(j)-V(50%)) are calculated, and furthermore a product 
of those difference values (xQ-l>V(50%))x ( x(j)-V(50%)) is calculated. 
When the x(t) crosses 50% level, i.e., zero-crossing level, the sign of the 
difference between the sample value and the V(50%), i.e., (x(j-l) -V(50%)) 
or (x(j)-V(50%)) changes from a negative sign to a positive sign or from a 

10 positive sign to a negative sign. Therefore, when the product is negative, it 
is detected that the x(t) has passed the zero-crossing level, and a time point j-1 
or j, at which a smaller absolute value of the difference between the sample 
value and the V (50%), i.e., (xG-1) -V(50%)) or (x(j)-V(50%)) is detected, is 
obtained as an approximated zero crossing point. Fig. 20 shows a waveform 

15 of a real part x(t) of an analytic signal. The marks o in Fig. 20 indicate the 
detected points (approximated zero-crossing points) closest to the 
corresponding rising zero-crossing points. 

WAVEFORM CLIPPING 

20 Waveform clipping means removes AM components from an input 

signal and retains only PM components in the input signal. A waveform 
clipping is performed by: 1) multiplying an analog or digital input signal by a 
constant, 2) replacing a signal value greater than a predetermined threshold 
value Thl with the threshold value Thl, and 3) replacing a signal value less 

25 than a predetermined threshold value Th2 with the threshold value Th2. 

Here, it is assumed that the threshold value Thl is greater than the threshold 
value Th2. Fig. 21 shows a clock signal having AM components. Since an 



-24- 

envelope of a time waveform fluctuates, it is seen that the AM components 
exist. Fig. 22 shows a clock signal that has been clipped by clipping means. 
Since the time waveform shows a constant envelope, it is confirmed that AM 
components have bee removed. 

5 

BRTEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 is a diagram showing an example of jitter measurement using a 
time interval analyzer; 

Fig. 2 is a diagram showing an example of histogram of period jitter 
1 0 measured using a time interval analyzer; 

Fig. 3 is a diagram showing an example of jitter measurement using a 

oscilloscope; 

Fig. 4 is a block diagram showing a functional configuration example 
of an interpolation based jitter measurement method; 
1 5 Fig. 5 A is a diagram showing a waveform example of a signal under 

measurement; 

Fig. 5B is a diagram showing a waveform example of period jitter of 
the signal under measurement shown in Fig. 5 A; 

Fig. 6 is a block diagram showing a functional configuration example 
20 of a j itter measurement by a conventional A§ method; 

Fig. 7A is a diagram showing a histogram example of a jitter 
measurement result by a time interval analyzer; 

Fig. 7B is a diagram showing a histogram example of a jitter 
measurement result by an interpolation based jitter measurement method; 
25 Fig. 8 is a diagram typically showing a calculation method of an 

instantaneous angular frequency; 

Fig. 9 is a diagram showing result examples of computer simulation 
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for respective period jitter estimations according to the method of the present 
invention in the case where a correction term is not used and in the case 
where a correction term is used; 

Fig. 1 OA is a diagram showing respective relationship examples 
5 between root-mean-square values of period jitter obtained by the method of 
the present invention and the numbers of over-samplings in the case of using 
the correction term and in the case of not using the correction term; 

Fig. 1 OB is a diagram showing respective relationship examples 
between peak-to-peak values of period jitter obtained by the method of the 
10 present invention and the numbers of over-samplings in the case of using the 
correction term and in the case of not using the correction term; 

Fig. 1 1 is a diagram showing an example of period jitter waveform of 
a clock signal under measurement; 

Fig. 12A is a diagram showing an example of period jitter histogram 
1 5 measured by the jitter measurement method according to the present 
invention; 

Fig. 12B is a diagram showing an example of period jitter histogram 
measured by a time interval analyzer; 

Fig. 13 is a diagram for comparing period jitter values measured by 
20 the time interval analyzer method with period jitter values measured by the 
jitter measurement method according to the present invention; 

Fig. 14 is a diagram showing an example of cycle-to-cycle period 
jitter waveform of a clock signal under measurement; 

Fig. 15 is a diagram showing an example of cycle-to-cycle period 
25 jitter histogram measured by the jitter measurement method according to the 
present invention; 

Fig. 1 6 is a diagram showing an example of digitized signal under 
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measurement; 

Fig. 17 is a diagram showing an example of power spectrum of a 
signal under measurement obtained by FFT; 

Fig. 1 8 is a diagram showing an example of band-limited power 
5 spectrum; 

Fig. 19 is a diagram showing an example of band-limited signal under 
measurement obtained by inverse FFT; 

Fig. 20 is a diagram showing an example of approximated 
zero-crossing points of a signal under measurement; 
10 Fig. 21 is a diagram showing an example of clock signal under 

measurement that has AM components; 

Fig. 22 is a diagram showing an example of clock signal under 
measurement that does not have AM components; 

Fig. 23 is a diagram showing an example of functional configuration 
15 of a jitter measurement apparatus according to the present invention; 

Fig. 24 is a flow-chart showing an example of jitter measurement 
method according to the present invention; 

Fig. 25 is a diagram showing another example of functional 
configuration of a jitter measurement apparatus according to the present 
20 invention; 

Fig. 26 is a flow-chart showing another example of jitter measurement 
method according to the present invention; 

Fig. 27 is a diagram showing an example of the functional 
configuration of the phase error estimator 104 in Fig. 23; 
25 Fig. 28 is a flow-chart showing an example of the phase error 

estimation method in the step 203 in Fig. 24; 

Fig. 29 is a diagram showing another example of the functional 
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configuration of the phase error estimator 104 in Fig. 23; 

Fig. 30 is a flow-chart showing another example of the phase error 
estimation method in the step 203 in Fig. 24; 

Fig. 3 1 is a diagram showing an example of the functional 
5 configuration of the period jitter estimator 105 in Fig. 23; 

Fig. 32 is a flow-chart showing an example of the period jitter 
estimation method in the step 204 in Fig. 24; 

Fig. 33 is a diagram showing another example of functional 
configuration of the period jitter estimator 105 in Fig. 23; 
10 Fig. 34 is a flow-chart showing another example of the period jitter 

estimation method the step 204 in Fig. 24; 

Fig. 35 is a diagram showing an example of the functional 
configuration of the band-pass filter 103 in Fig. 23; 

Fig. 36 is a flow-chart showing an example of the band-pass filtering 
1 5 method in the step 202 in Fig. 24; 

Fig. 37 is a diagram showing another example of the functional 
configuration of the band-pass filter 103 in Fig. 23; 

Fig. 38 is a flow-chart showing another example of the band-pass 
filtering method in the step 202 in Fig. 24; 
20 Fig. 39A is a diagram showing another configuration example of the 

input portion of a signal under measurement of a jitter measurement apparatus 
according to the present invention; 

Fig. 39B is a diagram showing further another configuration example 
of the input portion; 

25 Fig. 39C is a diagram showing further another configuration example 

of the input portion; 

Fig. 40A is a flow-chart showing a processing procedure of the input 
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portion shown in Fig. 39A; 

Fig. 40B is a flow-chart showing a processing procedure of the input 
portion shown in Fig. 39B; 

Fig. 41 is a diagram showing further another example of the period 
5 jitter estimator 105 in Fig. 23; and 

Fig. 42 is a flow-chart showing further another example of the period 
jitter estimation method in the step 204 in Fig. 24. 



DESCRIPTION OF THE EMBODIMENTS 

10 The embodiments of the present invention will be described below. 

Fig. 23 shows an example of configuration of a jitter measurement 
apparatus used in the embodiment of the present invention. This jitter 
measurement apparatus 100 comprises: a (fundamental) period estimator 101 
for obtaining a period T 0 of a fundamental wave of a signal under 

1 5 measurement; a memory 1 02 for storing therein the fundamental period 

estimated by the fundamental period estimator 101; a band-pass filter 103 for 
passing therethrough predetermined frequency components (components 
around the fundamental frequency) of the signal under measurement; a phase 
error estimator 104 for estimating sampling points (approximated 

20 zero-crossing points) t[k] (k=l ? 2, 3 5 — ) close to corresponding zero-crossing 
points of the signal that has passed through the band-pass filter 103, for 
obtaining phase errors 8[k] between the approximated zero-crossing points 
and the corresponding zero-crossing points of the signal under measurement, 
and for outputting a zero-crossing phase error data sequence 8[k] and a 

25 zero-crossing time interval sequence Ty^i between the approximated 

zero-crossing points; a period jitter estimator 105 for obtaining a period jitter 
sequence J[k] of the signal under measurement using the equation (26) from 
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the zero-crossing phase error data sequence 8[k] and the zero-crossing time 
interval sequence Tk,k+i estimated by the zero-crossing phase error estimator 
104 and the fundamental period To stored in the memory 102; and a jitter 
detector 106 to which the period jitter sequence J[k] is inputted for obtaining 
5 a jitter of the signal under measurement from the period jitter sequence J[k]. 

In addition, the jitter estimator 106 comprises a peak-to-peak detector 
107 for obtaining a difference between the maximum value and the minimum 
value of the period jitter sequence J[k] using the equation (29), an RMS 
detector 108 for calculating a root-mean-square value of the period jitter 

1 0 sequence J[k] using the equation (28), and a histogram estimator 109 for 
obtaining a histogram of the period jitter sequence J[k]. In addition, the 
band-pass filter 103 may be an analog filter or a digital filter, and moreover, 
the band-pass filter 1 03 may be a band-pass filter constituted by a software 
using FFT or the like as shown in Fig. 35 or Fig. 37. In addition, the phase 

15 error estimator 104 may be configured as shown in Fig. 27 or Fig. 29. In 
addition, the period jitter estimator 105 may be configured as shown in Fig. 
3 1 or Fig. 33 . A predetermined value can be stored in the memory 1 02 in 
advance as the fundamental period T 0 . In this case, the fundamental period 
estimator 101 can be omitted. 

20 Next, the operation in the case where a jitter measurement of the 

signal under measurement is performed using the jitter measurement 
apparatus 100 according to the present invention will be explained below. 
Fig. 24 shows a processing procedure of a jitter measurement method 
according to the present invention. First, in step 201, the fundamental period 

25 estimator 101 obtains the fundamental period T 0 of the signal under 

measurement, and the fundamental period T 0 is stored in the memory 102. 
Next, in step 202, the band-pass filter 103 selectively passes therethrough the 



-30- 

predetermined frequency components of the signal under measurement, i.e., 
frequency components around the fundamental frequency of the signal under 
measurement to band-limit the signal under measurement. Then, in step 
203, the phase error estimator 104 estimates the sampling points 
5 (approximated zero-crossing points) t[k] (k=l 9 2, 3, — ) close to the 

corresponding zero-crossing points of the signal that has passed through the 
band-pass filter 103, and obtains the phase errors 5[k] between the 
approximated zero-crossing points and the corresponding zero-crossing points 
of the signal under measurement to obtain the zero-crossing phase error data 

10 sequence 5[k] and the zero-crossing time interval sequence T kjk+ i between the 
corresponding approximated zero-crossing points. Next, in step 204, the 
period jitter estimator 105 estimates a period jitter sequence J[k] using the 
equation (26) from the zero-crossing phase error data sequence 8[k] and the 
zero-crossing time interval sequence T^k+i estimated by the zero-crossing 

1 5 phase error estimator 104 and the fundamental period T 0 stored in the memory 
102. Finally, in step 205, the jitter detector 106 obtains the period jitter of 
the signal under measurement from the period jitter sequence J[k], and the 
process ends. 

When the fundamental period T 0 is stored in the memory in advance, 
20 the step 201 for obtaining the fundamental period T 0 of the signal under 

measurement can be omitted. In addition, a processing procedure shown in 
Fig. 36 or Fig. 38 may be used as the step 202 for band-limiting the signal 
under measurement. In addition, a processing procedure shown in Fig. 28 or 
Fig. 30 can be used as the step 203 for estimating the zero-crossing phase 
25 error data sequence. In addition, a processing procedure shown in Fig. 32 or 
Fig. 34 may be used as the step 204 for estimating the period jitter of the 
signal under measurement. In addition, in the step 205 for obtaining the 
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period jitter of the signal under measurement, the peak-to-peak detector 107 
obtains the peak-to-peak value of period jitter using the equation (29), the 
RMS detector 108 obtains the RMS value of period jitter using the equation 
(28), and the histogram estimator 109 obtains the histogram from the 
5 instantaneous period waveform data. In order to estimate the approximated 
zero-crossing points, as apparent from the explanation of "approximated 
zero-crossing point detection method" in the aforementioned section of 
summary of the invention, the signal to which a band-pass filtering process 
has been applied must be an over-sampled (sampling at a frequency equal to 

10 or higher than three times of the fundamental frequency) digital signal. If 
the output of the band-pass filter 103, i.e., the band-limited signal in the step 
202 is an analog signal, the signal inputted to the phase error estimator 104 is 
converted into a digital signal within the phase error estimator 104, whereby 
the digital signal is processed. That is, in the step 203, the signal is 

15 processed after the input signal is converted into a digital signal. 

Fig. 25 shows another example of configuration of a jitter 
measurement apparatus used in the embodiment of the present invention. 
This jitter measurement apparatus 300 is same as the jitter measurement 
apparatus shown in Fig. 23 except that the jitter measurement apparatus 300 

20 includes a cycle-to-cycle period jitter estimator 301 to which the period jitter 
sequence S[k] is inputted for calculating a difference waveform using the 
equation (30) to output a cycle-to-cycle period jitter sequence JccM, and a 
switch 302 for selecting an output sequence from the period jitter estimator 
106 or the cycle-to-cycle period jitter estimator 301 . For simplicity, the 

25 explanation of the overlapped portion is omitted. In this case, the jitter 
detector 106 obtains a jitter value (a peak-to-peak value, an RMS value or a 
histogram) of the signal under measurement from the period jitter sequence 
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J[k] or the cycle-to-cycle period jitter sequence JccM- 

Next, the operation in the case where a jitter measurement of the 
signal under measurement is performed using the jitter measurement 
apparatus 300 according to the present invention will be explained below. 
5 Fig. 26 shows another processing procedure of a jitter measurement method 
according to the present invention. This jitter measurement method is same 
as the jitter measurement method shown in Fig. 24 except that this jitter 
measurement method includes: a step 401 for calculating, from the period 
jitter sequence J[k] estimated by the period jitter estimator 105, a difference 

10 waveform of the period jitter sequence J[k] by means of the cycle-to-cycle 
period jitter estimator 301, after the period jitter of the signal under 
measurement is obtained from the period jitter sequence J[k], and for 
obtaining a cycle-to-cycle period jitter sequence JccM; and a step 402 for 
obtaining a cycle-to-cycle period jitter of the signal under measurement from 

15 the cycle-to-cycle period jitter sequence Jcc[k] in the state that the jitter 
detector 106 connects the switch 302 to the cycle-to-cycle period jitter 
estimator 301 side. In this case, the explanation of the overlapped portion is 
omitted. In the step 402 for obtaining the cycle-to-cycle period jitter of the 
signal under measurement, the peak-to-peak detector 107 obtains a 

20 peak-to-peak value of period jitter using the equation (32), the RMS detector 
108 obtains an RMS value of period jitter using the equation (31), and the 
histogram estimator 109 obtains a histogram from the instantaneous period 
sequence. 

The jitter measurement apparatus shown in Fig. 25 can also be 
25 constructed as an apparatus for estimating only cycle-to-cycle period jitter. 
In this case, the switch 302 for selecting a sequence can be omitted. 
Similarly, the jitter measurement method shown in Fig. 26 may also estimate 
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only cycle-to-cycle period jitter. In this case, the step 205 for estimating the 
period jitter value from the period jitter sequence 5[k] can be omitted. 

Fig. 27 shows an example of configuration of the phase error 
estimator 104 used in the jitter measurement apparatus 100. The phase error 
5 estimator 104 comprises: a zero-crossing sampler 501 to which a signal that 
has passed through the band-pass filter 103 is inputted for sampling only 
waveform data close to zero-crossing timings of the band-limited signal, and 
for outputting an approximated zero-crossing data sequence Xzx[k]=x(t[k]) 
and a zero-crossing time interval sequence T k)k+ i between the approximated 

1 0 zero-crossing points t[k] and t[k+l]; and a phase error calculator 502 to which 
the approximated zero-crossing data sequence x zx [k] is inputted for 
calculating phase errors 8[k] between the approximated zero-crossing points 
t[k] and the corresponding zero-crossing points of the signal under 
measurement using the equation (18) or (20) from the approximated 

15 zero-crossing data x Z x[k] ? and for outputting a zero-crossing phase error data 
sequence 5[k]. 

Next, Fig. 28 shows a processing procedure for estimating the phase 
errors between the zero-crossing points of the signal under measurement and 
the corresponding approximated zero-crossing points using the phase error 

20 estimator 104 according to the present invention. First, in step 601, the 
zero-crossing sampler 501 samples only the waveform data close to the 
zero-crossing timings of the signal that has passed through the band-pass filter 
103 to obtain the approximated zero-crossing data sequence x zx [k]==x(t[k]). 
Next, in step 602, the zero-crossing sampler 501 obtains the approximated 

25 zero-crossing timing sequence t[k] of zero-cross-sampled approximated 

zero-crossing data, and calculates a difference sequence of the approximated 
zero-crossing timing sequence t[k] to obtain the zero-crossing time interval 
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sequence T k)k+ i=t[k+l]-t[k]. Finally, in step 603, the phase error calculator 
502 calculates the phase errors between the approximated zero-crossing 
points and the corresponding zero-crossing points of the signal under 
measurement from the approximated zero-crossing data sequence x zx [k] 
5 obtained by the zero-crossing sampler 50 1, and outputs the zero-crossing 
phase error data sequence 5[k], and then the process ends. In the step 602 
for obtaining the zero-crossing time interval sequence, the zero-crossing 
sampler 501 obtains the zero-crossing time interval T^k+i between the 
approximated zero-crossing points using the equation (11). In addition, in 

10 the step 603 for obtaining the zero-crossing phase error data sequence, the 
phase error calculator 502 obtains the phase errors between the approximated 
zero-crossing points and the corresponding zero-crossing points using the 
equation (1 8) or (20). In the case of obtaining the phase errors using the 
equation (1 8) or (20), if the signal under measurement x(t) is not a sine wave, 

15 it is necessary to make it a sine wave using the band-pass filter 103. Further, 
the amplitude A in the equations (18) and (20) is obtained by calculating a 
root-mean-square value of output samples of the band-pass filter 103 or by 
calculating one half of a difference between the maximum value and the 
minimum value of the output samples of the band-pass filter 103. 

20 Fig. 29 shows another example of configuration of the phase error 

estimator 104 used in the jitter measurement apparatus 100. This phase error 
estimator 104 is same as the phase error estimator 104 shown in Fig. 27 
except that this phase error estimator 104 includes an analytic signal 
transformer to which the signal that has passed through the band-pass filter 

25 103 is inputted for transforming the band-limited signal into a complex 
analytic signal. In this case, the explanation of the overlapped portion is 
omitted. In this case, the analytic signal is inputted to the zero-crossing 
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sampler 501, by which analytic signal waveform data around zero-crossing 
timings of a real part of the analytic signal are sampled, and an approximated 
zero-crossing complex data sequence z(t[k])=XzxM + j *zx[k] (equation (22)) 
is estimated. In addition, the phase error calculator 502 obtains a 
5 zero-crossing phase error data sequence 8[k] by calculating the equations (23) 
and (24) from the approximated zero-crossing complex data sequence. 

Next, a processing procedure for estimating the phase errors between 
the zero-crossing points of the signal under measurement and the 
corresponding approximated zero-crossing points using the phase error 

10 estimator 104 shown in Fig. 29 will be explained with reference to Fig. 30. 
First, in step 801, the analytic signal transformer 701 transforms the 
band-limited signal obtained from the band-pass filter 103 into a complex 
analytic signal z(t). Next, in step 802, the zero-crossing sampler 501 
samples the analytic signal waveform data close to the zero-crossing timings 

15 of the real part x(t) of the analytic signal z(t) transformed by the analytic 
signal transformer 701 to estimate the approximated zero-crossing complex 
data sequence z(t[k]). Next, in step 803, the zero-crossing sampler 501 
obtains a timing sequence t[k] of the sampled zero-crossing complex data, and 
calculates a difference sequence of the timing sequence t[k] to obtain the 

20 zero-crossing time interval sequence T k5k+ i=t[k+l]-t[k]. Finally, in step 804, 
the phase error calculator 502 calculates the phase errors between the 
approximated zero-crossing points and the corresponding zero-crossing points 
of the signal under measurement from the approximated zero-crossing 
complex data sequence z(t[k]) obtained by the zero-crossing sampler 501, and 

25 obtains the zero-crossing phase error data sequence 8[k], and then the process 
ends. 

In the step 801 for transforming the band-limited signal into a 
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complex analytic signal z(t), the analytic signal transformer 701 transforms an 
input signal into an analytic signal using Hilbert Transform or Fast Fourier 
Transform. In addition, in the step 803 for obtaining a zero-crossing time 
interval sequence, the zero-crossing sampler 501 obtains zero-crossing time 
5 intervals Tk,k+i between the approximated zero-crossing points using the 
equation (11). In addition, in the step 804 for obtaining a zero-crossing 
phase error data sequence, the phase error calculator 502 obtains the 
instantaneous phases <f>[k] of the approximated zero-crossing points using the 
equation (23), and obtains the phase errors 5[k] between the approximated 

10 zero-crossing points and the corresponding zero-crossing points using the 
equation (24) or (25). 

Fig. 3 1 shows an example of configuration of the period jitter 
estimator 105 used in the jitter measurement apparatus 100. This period 
jitter estimator 105 comprises: an instantaneous phase estimator 901 for 

1 5 obtaining an instantaneous period sequence T[k] of the signal under 

measurement from the zero-crossing phase error data sequence S[k], the 
zero-crossing time interval sequence T^k+u and the fundamental period To of 
the signal under measurement; a fundamental period subtracter 902 for 
obtaining a difference between each instantaneous period value T[k] 

20 estimated by the instantaneous period estimator 901 and the fundamental 

period To, and for outputting a period jitter sequence (uncorrected period jitter 
sequence); a period jitter sequence corrector 903 for correcting the period 
jitter sequence, which is estimated using difference formula, by multiplying 
the period jitter sequence (uncorrected period jitter sequence) by a ratio 

25 (T 0 /Tk,k+i) of the fundamental period T 0 and the zero-crossing time interval 
Tk,k+i- 

Next, the operation in the case of performing a period jitter estimation 
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of the signal under measurement using this period jitter estimator 105, i.e., a 
processing procedure of the period jitter estimation method will be explained 
with reference to Fig. 32. First, in step 1001, the instantaneous period 
estimator 901 obtains an instantaneous period sequence T[k] from the 
5 zero-crossing phase error data sequence 5[k] estimated by the zero-crossing 
phase error estimator 104, the zero-crossing time interval sequence T k)lc+ i 
obtained from the zero-crossing sampler 501, and the fundamental period T 0 
of the signal under measurement stored in the memory 102. Next, in step 
1002, the fundamental period subtractor 902 calculates a difference between 

1 0 each instantaneous period value T[k] estimated by the instantaneous period 
estimator 901 and the fundamental period To to obtain a period jitter 
sequence. Finally, in step 1003, the period jitter sequence corrector 903 
corrects the period jitter sequence, which is estimated using difference 
formula, by multiplying the period jitter sequence obtained by the 

15 fundamental period subtractor 902 by a ratio (T 0 /T k)k+ i) of the fundamental 
period and the zero-crossing time interval, and then the process ends. In the 
step 1001 for obtaining the instantaneous period sequence T[k] of the signal 
under measurement, the instantaneous period estimator 901 obtains the 
instantaneous period T[k] using the equation (27). 

20 Fig. 33 shows another example of configuration of the period jitter 

estimator 105 used in the jitter measurement apparatus 100. This period 
jitter estimator 105 comprises: an instantaneous angular frequency estimator 
1001 to which the zero-crossing phase error data sequence 8[k] is inputted for 
obtaining an instantaneous angular frequency sequence o[k] by a calculation 

25 of the equation (10) from the sampling intervals Tk,k+i of the approximated 
zero-crossing data sequence and the zero-crossing phase error data 8[k]; and a 
period jitter calculator 1 102 for obtaining a period jitter sequence J[k] by a 
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calculation of the equation (9) from the instantaneous angular frequency 
sequence to[k] and the fundamental period To. 

Next, the operation in the case of performing a period jitter estimation 
of the signal under measurement using this period jitter estimator 105 shown 
5 in Fig. 33 will be explained. Fig. 34 shows a processing procedure of this 
period jitter estimation method. First, in step 1201, the instantaneous 
angular frequency estimator 1101 obtains the instantaneous angular frequency 
sequence a>[k] of the signal under measurement from the zero-crossing phase 
error data 5[k] estimated by the zero-crossing phase error estimator 104 and 

10 the zero-crossing time interval sequence T kjk+1 obtained from the 

zero-crossing sampler 501. Next, in step 1202, the period jitter calculator 
1 102 obtains the period jitter sequence J[k] of the signal under measurement 
from the instantaneous angular frequency sequence co[k] obtained from the 
instantaneous angular frequency estimator 1101 and the fundamental period 

1 5 T 0 stored in the memory 102, and the process ends. In the step 1201 for 
obtaining the instantaneous angular frequency sequence co[k] of the signal 
under measurement, the instantaneous angular frequency estimator 1101 
obtains the instantaneous angular frequency using the equation (10). In 
addition, in the step 1202 for obtaining the period jitter J[k] of the signal 

20 under measurement from the instantaneous angular frequency sequence oi>[k], 
the period jitter calculator 1102 obtains the period jitter sequence J[k] using 
the equation (9). 

Fig. 35 shows an example of configuration of the band-pass filter 103 
used in the jitter measurement apparatus 100. This band-pass filter 103 
25 comprises: a time domain to frequency domain transformer 1301 for 

transforming the signal under measurement into a signal in frequency domain; 
a bandwidth controller 1302 for extracting only components around the 
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fundamental frequency of the signal under measurement from the output of 
the time domain to frequency domain transformer 1301; and a frequency 
domain to time domain transformer 1303 for inverse-transforming the output 
of the bandwidth controller 1302 into a signal in time domain. The time 
5 domain to frequency domain transformer 1 30 1 and the frequency domain to 
time domain transformer 1303 may be packaged using FFT and inverse FFT, 
respectively. 

Next, the operation in the case of performing a band-limitation of the 
signal under measurement using this band-pass filter 103 shown in Fig. 35 

1 0 will be explained. Fig. 36 shows a processing procedure of this 

band-limitation method. First, in step 1 40 1 , the time domain to frequency 
domain transformer 1301 applies FFT to the signal under measurement to 
transform a signal in time domain into a signal in frequency domain. Next, 
in step 1402, the bandwidth limiter 1302 retains the components around the 

1 5 fundamental frequency of the signal under measurement in the transformed 
signal in frequency domain, and replaces the other frequency components 
with zero to band-limit the signal in frequency domain. Finally, in step 
1403, the frequency domain to time domain transformer 1303 applies inverse 
FFT to the band-limited signal in frequency domain to transform the signal in 

20 frequency domain into a signal in time domain, and the process ends. 

Fig. 37 shows another example of configuration of the band-pass filter 
1 03 used in the jitter measurement apparatus 1 00. This band-pass filter 1 03 
comprises: a buffer memory 1501 for storing therein the signal under 
measurement; a data selector 1502 for extracting the signal in the sequential 

25 order from the buffer memory 1501 such that the signal being extracted is 
partially overlapped with the signal extracted just before; a window function 
multiplier 1503 for multiplying each extracted partial signal by a window 
function; a time domain to frequency domain transformer 1504 for 
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transforming each partial signal multiplied by the window function into a 
signal in frequency domain; a bandwidth limiter 1505 for extracting only 
components around the fundamental frequency of the signal from the signal 
transformed in frequency domain; a frequency domain to time domain 
5 transformer 1 506 for inverse-transforming the output of the bandwidth limiter 
1505 into a signal in time domain; and an amplitude corrector 1507 for 
multiplying the signal transformed in time domain by an inverse number of 
the window function to obtain a band-limited analytic signal. The time 
domain to frequency domain transformer 1504 and the frequency domain to 
1 0 time domain transformer 1 506 may be packaged using FFT and inverse FFT, 
respectively. 

Next, the operation in the case of performing the band-limitation of 
the signal under measurement using this band-pass filter 103 will be 
explained. Fig. 3 8 shows a processing procedure of this band-limitation 

1 5 method. First, in step 1601, the buffer memory 1501 stores therein the 
signal under measurement. Next, in step 1602, the data selector 1 502 
extracts a portion of the stored signal from the buffer memory 1501. Next, 
in step 1603, the window function multiplier 1503 multiplies the extracted 
partial signal by the window function. Next, in step 1604, the time domain 

20 to frequency domain transformer 1 504 applies FFT to the partial signal 

multiplied by the window function to transform the signal in time domain into 
a signal in frequency domain. Next, in step 1 605, the bandwidth limiter 
1505 retains only components around the fundamental frequency of the signal 
under measurement in the transformed signal in frequency domain, and 

25 replaces the other frequency components with zeros to band-limit the signal in 
frequency domain. Then in step 1 606, the frequency domain to time domain 
transformer 1506 applies inverse FFT to the band-limited signal in frequency 
domain to transform the signal in frequency domain into a signal in time 
domain. Next, in step 1607, the amplitude corrector 1507 multiplies the 
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inverse-transformed signal in time domain by an inverse number of the 
window function multiplied in the step 1603 to obtain a band-limited signal 
Finally, in step 1608, a check is made to see if there is unprocessed data 
remaining in the buffer memory. If there is unprocessed data remaining in 
5 the buffer memory, the data selector 1502 extracts, in step 1609, the signal in 
the sequential order from the buffer memory 1501 such that the signal being 
extracted is partially overlapped with the signal extracted just before, and 
thereafter the steps 1603, 1604, 1605, 1606, 1607 and 1608 are repeated. If 
there is no unprocessed data, the process ends. 

10 The band-pass filter 103 shown in Fig. 37 and the band-pass filtering 

process shown in Fig. 38 are modified versions of the band-pass filter 103 
shown in Fig. 35 and the band-pass filtering process shown in Fig. 36, 
respectively- This modification is made to perform a real time processing of 
the band-pass filter 103 shown in Fig. 35 and the band-pass filtering process 

15 shown in Fig. 36. Each of the band-pass filters 103 shown in Figs. 35 and 
37 commonly has a time domain to frequency domain transformer, a 
bandwidth limiter and a frequency domain to time domain transformer, and 
each of the band-pass filtering processes shown in Figs. 36 and 38 commonly 
includes a time domain to frequency domain transformation step, a bandwidth 

20 limitation step and a frequency domain to time domain transformation step. 

In the jitter measurement apparatus shown in Figs. 23 and 25, as 
shown in Fig. 3 9 A, an analog signal under measurement may be digitized and 
converted to a digital signal by the AD converter 1701, and then the digitized 
signal may be supplied to the period estimator 101 and the band-pass filter 

25 103. It is desirable that a high-speed AD converter, a high-speed digitizer or 
a high-speed digital sampling oscilloscope is used as the AD converter 1701. 
Regarding the jitter measurement method corresponding to this, prior 
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to the estimation of the fundamental period T 0 in the step 201 in Figs. 24 and 
26, it is sufficient to sample (digitize) the analog signal under measurement 
and to convert the sampled signal to a digital signal in step 1801 as shown in 
Fig. 40A. 

5 In the jitter measurement apparatus shown in Figs. 23 and 25, as 

shown in Fig. 39B, amplitude modulation components may be removed from 
the analog signal under test by the waveform clipper 1901 so that the signal is 
in the state that the phase modulation components are retained, and then the 
signal may be supplied to the period estimator 101 and the band-pass filter 
10 103. 

Regarding the jitter measurement method corresponding to this, prior 
to the estimation of the fundamental period T 0 in the step 201 in Figs. 24 and 
26, it is sufficient to remove, in step 2001 as shown in Fig.40B, the amplitude 
modulation components from the input signal so that it is in the state that the 

15 phase modulation components are retained. 

As shown in step 39C, the AD converter 1701 may be provided at the 
rear stage or at the front stage of the waveform clipper 1901 so that the 
amplitude modulation components of the analog signal under measurement 
are removed from the signal, and the analog signal under measurement is 

20 converted into a digital signal, and then the digital signal is supplied to the 
period estimator 101 and the band-pass-filter 103. Although not shown, in 
the jitter measurement method according to the present invention, the 
fundamental period T 0 of the signal under measurement may be similarly 
estimated after the amplitude modulation components of the signal under 

25 measurement are removed and the signal under measurement is converted into 
a digital signal. In this case, either of the removing process of the amplitude 
modulation components and the analog to digital conversion process may be 
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performed first. 

Fig. 41 shows further another example of configuration of the period 
jitter estimator 1 05 used in the embodiment of the present invention. This 
period jitter estimator 105 is same as the period jitter estimator 105 shown in 

5 Fig. 3 1 except that this period jitter estimator 105 includes a low-pass filter 
2101 for removing high frequency components from the estimated period 
j itter sequence outputted by the period j itter sequence corrector 903 . In this 
case, the explanation of the overlapped portion is omitted. 

Fig. 42 shows a processing procedure of a period jitter estimation 

1 0 method corresponding to the period jitter sequence estimator 1 05 shown in 
Fig. 4 1 . This period jitter estimation method is same as the period jitter 
estimation method shown in Fig. 32 except that this period jitter estimation 
method includes, in the final step, a step 2201 in which the low-pass filter 
2101 removes high frequency components from the period jitter sequence. 

1 5 As shown in Figs. 41 and 42, the removal of the high frequency 

components by the low-pass filter 2101 is intended to remove the components 
close to a high-frequency noise component like the fundamental frequency 
1/T 0 of the signal under measurement, and hence the cut-off frequency is, for 
example, 1/2 of the fundamental frequency, i.e., 1/(2T 0 ) or so. 

20 The low-pass filter 2 1 0 1 can also be built, as indicated by dashed lines 

in Fig. 3 3 , in the period jitter estimator 105. In this case, the processing 
procedure of the period jitter estimation method shown in Fig. 34 includes, in 
the final stage of the process, a step of removing high-frequency components 
of a period jitter sequence using the low-pass filter. 

25 In the above description, if the signal under measurement is a signal 

close to a sine wave, the band-pass filter 103 and the band-pass filtering 
process 202 may be omitted. 
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In Figs. 23 and 25, the jitter detector 106 may include at least one of 
the peak-to-peak detector 107, the RMS detector 108 and the histogram 
estimator 1 09. Similarly, in Figs. 24 and 26, each of the period jitter value 
estimation process and the cycle-to-cycle period jitter value estimation 
5 process may include at least one of the peak-to-peak detection process, the 
root-mean-square process and the histogram estimation process. 

The jitter measurement apparatus according to the present invention 
shown in Figs. 23 and 24 may be functioned by the execution of programs in 
a computer. 

1 0 According to the j itter measurement apparatus and the j itter 

measurement method according to the present invention, the estimation 
accuracy in the case where the signal under measurement is digitized at low 
over-sampling rate can be improved. Therefore, the accuracy of jitter 
measurement using the conventional A<j) method can greatly be improved. 



